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Abstract 



Systems of partial differential equations lie at the heart of physics. De- 
spite this, the general theory of these systems has remained rather obscure 
in comparison to numerical approaches such as finite element models and 
various other discretisation schemes. There are, however, several theoret- 
ical approaches to systems of PDEs, including schemes based on differ- 
ential algebra and geometric approaches including the theory of exterior 
differential systems j|] and the so-called "formal theory" Q built on the 
jet bundle formalism. This paper is a brief introduction to jet bundles 
focusing on the completion of systems to equivalent involutive systems for 
which power series solutions may be constructed order by order. We will 
not consider the mathematical underpinnings of involution (which lie in 
the theory of combinatorial decompositions of polynomial modules [|) 
nor other applications of the theory of jet bundles such as the theory 
of symmetries of systems of PDEs pi or discretisation schemes based on 
discrete approximations to jet bundles jj]] . 

1 Fibre Bundles and Sections 

A bundle is a triple (M, X, tt), where M is a manifold called the total space, X 
is a manifold called the base space, and tt : M — > X is a continuous surjective 
mapping called the projection. Where no confusion can arise, we shall often find 
it convenient to denote the bundle either by its total space or its projection. A 
trivial bundle is a bundle whose total space M is homeomorphic to X X U, where 
U is a manifold called the fibre. A bundle which is locally a trivial bundle is 
called a fibre bundle (an example of a fibre bundle which is not trivial is the 
Mobius band). In the following we shall only be concerned with trivial bundles. 

We will denote the coordinates on the base manifold x — {x l , i — 1, . . . ,p} 
where p is the dimension of the base manifold and the coordinates on the fibre 
u = {u a , a = 1, . . . , q} where q is the dimension of the fibre (when we consider 
jet bundles, we will need to extend this notation slightly). In the case of a trivial 
bundle, or in a local coordinate patch on a fibre bundle, the projection takes 
the simple form 




A section of a fibre bundle is a map 
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such that 7r o <j>j is the identity map on X. In other words, a section assigns to 
each point in X a point in the fibre over that point. The graph of the function 
f(x) is: 

T f = {(x,f(x)) :xen}dXxU 

where fi is the domain of definition of /. We will find it convenient to refer to 
sections, functions and graphs interchangably. 

2 Jet Bundles 

We define a multi-index J as a p-tuple [ji, J2, • ■ • ,j p ] with ji G N. The order of 
the multi- index J, denoted \J\ is given by the sum of the jj. We will often find 
it more convenient to use a repeated-index notation for J. In this notation J is 
represented by a string of \J\ independent coordinate labels, with ji copies of 
the z-th coordinate label. For example, if p — 3 and the coordinates are labelled 
then x, y and z then the second order multi-indices in repeated index notation 
are xx, xy, xz, yy, yz and zz. We introduce the special notation J, i where i is an 
independent coordinate label for the multi-index given by . . . , ji + 1, . . . ,j p }. 
For example, xyy, x = xxyy. 

If our independent variables are x % and our dependent variables are u a then 
we introduce jet variables Uj where J is a multi-index. Notice that we can put 
the jet variables of order n in one-to-one correspondance with the derivatives 
of the dependent variables of order n. We will later introduce further struc- 
tures that enable us to make a full correspondance between jet variables and 
derivatives. 

Associated with these jet variables we introduce a set of Euclidean spaces 
Ui, whose coordinates are Uj with \J\ — i. We call the space M« = X x U x Ui 
the first order jet bundle over the space M = X x U. We now introduce the 
notation 

U {n) =U Xf/lX • • • xU n 

and call the space = X x f/M the n-th order jet bundle over M. 

Example 1 Let p = 2 and q = 1. Label the independent variables x and 
y and the dependent variable u. The first order jet bundle, then has 

coordinates (x,y,u,u x ,u y ), the second order jet bundle, M^ 2 ' , has coordinates 
(x,y,u,u x ,u y ,u xx ,u xy ,Uyy) and so on. 

We will often consider a jet bundle as a bundle over a lower order jet bundle. 
We denote the natural projection between the (m + n)-th order jet bundle and 
the n-th order jet bundle as 

7r™ +Tl : M^ m+ ™' — > 

Note that although M^ m+n ^ is a bundle over it is not a jet bundle over 

M^ n \ but rather a subset of such a bundle. 
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3 Differential Functions and Formal Derivatives 



A differential function is a smooth, real- valued function defined on M^ n \ We 
denote the algebra of differential functions defined on M^> by A( n y If F G A( n ) 
then F e A( m+n ) too, as the coordinates on are a subset of the coordinates 
on M( m+n \ If the lowest order space on which F is defined is then we will 
say that F is an n-th order differential function. These will be used to describe 
sections of M and differential equations. 

The most fundamental maps between lower order and higher order jet bun- 
dles are provided by formal derivatives. A formal derivative operator Di, called 
the formal derivative with respect to x l , maps each differential function F e A( n ) 
to a differential function DiF e -4(„+i) via 

a = l J •' 

It is convenient to extend the notation for formal and partial derivatives to 
encompass our multi-index notation: 

D J = (D 1 )*(W 2 ...(A,) A ' 
and similarly for partial derivatives. Clearly if F £ A( n ) then DjF G A( n +\j\y 



4 Prolongation of Sections 

If Tf cMisa section defined by u = f(x) where / is a smooth function of x 
then we can use the formal derivative to prolong it to a section T^ 1 ' C M^. 

The equations defining are simply u a = f a (x) and the equations found by 
applying each of the Di to u a = f a (x): 




= rw 

= d t f a (x) 



Similarly we can use the Di multiple times to prolong Tf to a section 1^ C 
M( n ) defined by the equations 

U a j = djf a {x) 

where J ranges over all multi-indices such that ^ \ J\ ^ n. We will sometimes 
talk about the n-th prolongation of a function f(x), and write this prolongation 

as f {n) (x). 



5 Differential Equations and Solutions 

We intend to view systems of PDEs as geometric objects. In keeping with this 
programme, we will simply call such a system a "differential equation" . An n-th 
order differential equation, I a, is a fibred submanifold of M^ n \ The differential 
equation is often stated as the kernel of a set of differential functions A v £ A( n ) : 

A„(x,u^) = 0, v =l,...,l 
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We can map any system of partial differential equations onto such a subman- 
ifold, simply by replacing all of the derivatives of dependent variables by the 
corresponding jet variables. Indeed we have chosen our notation in such a way 
that this process is entirely transparent. 

Example 2 The two dimensional wave equation 

d 2 u d 2 u d 2 u 

dt 2 dx 2 dip- 
maps onto the submanifold of (with the obvious coordinates) determined 
by 

Utt l^xx ^yy — 
In keeping with our compact notations x, u and /, we shall write 

A = (A!(s, «<">),..., A,(s, «<">)) eA\ n) 

A is therefore a map from to R ( . The differential equation 7a is the 

submanifold in which the map A vanishes: 

I A = {(i, U W) : A(i,uW) = 0} 

A smooth solution of 7a C is a smooth function f{x) such that 

A(x,/(")(x))=0 
or, in terms of our geometric formulation 

Tf c 7 A 

Note that not every section of M^ n > which lies entirely within 7a is a prolon- 
gation of a section of M - the prolongation of a section of M automatically 
respects the correspondence between jet variables and derivatives, whereas an 
arbitrary section of M^> does not. 

6 Prolongation and Projection of Differential 
Equations 

The fc-th prolongation of the differential equation 

7 A = {(x,u^) : A„(jc,u( n >) = 0} C AfW 

is 

7^ fe) = {(x, U (" +fe )) : (Dj/\ v )(x,u {n+k) ) = 0} C M(" +fe ) 

where J runs over all multi-indices up to order k. 

Differential equations may be projected along the fibers onto lower order jet 
bundles. In general, this is a complicated procedure in local coordinates, but 
it is much easier if the differential equation is known to be the prolongation of 
a lower order system. For the remainder of this paper this will always be the 
case. 
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Example 3 Let the differential equation I C be defined by 



U ZZ + U X y + U = 

Ux — u = 

Uy — U 2 = 



The projection of 2 into is 

7T?(Z) 



u = 
,i 2 = 



7 Power Series Solutions 

A smooth solution of I C Af( n ) in a neighbourhood of xo may be written as the 
power series 

« Q = / Q (ao = £ £{x- Xo y 

\J\=0 

for some constants a". Here J! = ji! J2! • ■ • j p ! and 

i=l 

All we have to do is choose the values of the a" so that C X. It will, however, 

prove to be easiest to fix the a" by using the condition that r^" +fe) C X {k) for 
all k and working entirely at the point x - By applying the formal derivative 
repeatedly to the power series we can obtain power series expressions for each 
of the Uj(x). Evaluating each of these at xo shows that a" — Uj\ Xg . Therefore 
we require that 

DjA v (xo,a%)=0 

for all J. We have exchanged the solution of a set of partial differential equations 
for the solution of an infinite number of algebraic equations. 

For a class of systems known as formally integrable differential equations we 
can construct a power series solution order by order. We first substitute the 
general form for the power series into each of the equations and evaluate at 
x = xq- This gives us a set of algebraic equations for the a"(x) with \ J\ ^ n. 
We then make a partition of the jet variables into parametric derivatives whose 
values we can choose and principal derivatives whose values are then fixed by 
the system, and solve for the latter in terms of the former. We then prolong 
the differential equation and repeat the process. This time the equations of 
order less than n + 1 will automatically be satisfied by the previously chosen 
constants, and we will be left with a new set of equations DiA l/ (xo, a") = for 
I J I = n + 1. The nature of the formal derivative means that these equations 
will be linear. We may repeat the procedure to calculate ever higher terms in 
the power series. 
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8 Integrability Conditions 

If a differential equation is not formally integrable the solutions are subject to 
constraints which we call integrability conditions. These are extra equations 
that are differential rather than algebraic consequences of the equations A = 0. 
In other words, projecting the prolongation of a differential equation may not 
return the original equation but only a proper subset thereof: 

and so the order by order construction of a power series solution will be dis- 
rupted. To streamline the notation, we will write the j-th projection of the fc-th 
prolongation of 2 as For example, the expression above may be rewritten 

4 fe) £ I- 

Integrability conditions arise in two ways: through the differentiation of 
equations with order less than ninlc M^ n \ and through the effects of cross- 
derivatives, as shown in the following example: 

Example 4 Let p = 3 with coordinates x, y and z, and q = 1 with coordinate 
u, and consider the differential equation 

2 ( u z +yu x = 
• \ u y = 

which prolongs to 

u xz + yu xx = 
u yz + u x + yu xy = 
u zz + yu xz = 

u yy = 
u yz = 
u z + yu x = Q 
u y = 

We see that the equations u yz = and u xy = substituted into the equation 
u yz + u x + yu xy = imply that u x = 0. This is a first order equation and so 
forms part of on projection. Hence £ 1. 

A differential equation X C M^> is formally integrable if x[ k+1 ^ = for 
all k. Notice that to check for formal integrability requires an infinite number of 
operations, for integrability conditions may in general arise after an arbitrarily 
large number of prolongations. 



T« : <^ 



9 Involutive Differential Equations 

We now turn to the consideration of a subset of formally integrable differential 
equations known as involutive equations. Two facts make this class of equations 
interesting and useful. Firstly, it is possible to determine whether a given differ- 
ential equation is involutive using only a finite number of operations. Secondly, 
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for any differential equation it is possible to produce an involutive equation with 
the same solution space using only a finite number of operations. 

There is a more systematic method for determining the intcgrability condi- 
tions that arise upon a single prolongation and projection. Let us look at the 
Jacobi matrix of X^ C M^ n+1 \ This matrix can be divided into four blocks: 





W ' |J|_n + 1 








We order the columns according to increasing | J|, and within each order by the 
first non- vanishing component of J (we will call this component the class of J). 

When we project the X^ into M (") to form jW we must include only those 
equations that are independent of the Uj with J| = n + 1. In other words 
we must include all those equations which have a full row of zeros in the right 
hand block. This clearly includes the equations corresponding to rows in the 
bottom part of the Jacobi matrix. However, if the upper right submatrix is not 
of maximal rank then we may be able to form intcgrability conditions. If for 
a row with all zeros in the right hand section we find that the left hand part 
is independent of the rows in the lower part of the matrix then there is indeed 
an intcgrability condition, which can be determined by performing the same 
operations on the full equations DiA v = 0. 

We will call the system of equations defined by the upper right block of the 
Jacobi matrix the symbol of X, and denote this Sym X: 



where the Vj are a new set of variables, which we order in the same way as 
the the Uj when displaying the symbol as a matrix. Notice that the entries in 
the matrix of Sym 2" are the coefficients of the highest order jet variables in the 
equations defining X^ as can be seen by comparison with the formal derivative. 

Comparison of the ranks of X, X^ and SymZ^ 1 ) will enable us to determine 
if an intcgrability condition will occur on a single prolongation and projec- 
tion. There is an intcgrability condition if ran rankZ, or equivalently if 
dimZj 1 ^ < dimZ. Furthermore from inspection of the Jacobi matrix of X^> 



We can thus systematically determine if integrability conditions arise from a 
single prolongation, and if necessary find the new equations. 

Henceforth we will always consider the row echelon form of the symbol. We 
call Sym X involutive if 



where (3k is the number of rows of class k in Sym X. 

For a row of class k we call the variables x 1 , x 2 , . . . , x k multiplicative vari- 
ables. We now consider prolonging each equation by its multiplicative variables 




ranklf ) = rankT (1) - rankSym X {1) 




k 
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only. The equations obtained in this manner will be independent as they will 
have distinct pivots in Sym X^ . As there are (3k equations of class k and each 
has k multiplicative variables then this means there will be at least ^ kf3k in- 
dependent equations of order n + 1 in l( n+1 \ If Sym I is involutive then we 
obtain all the independent equations of order n + 1 in this manner. The equa- 
tions obtained from the other prolongations required to prolong X to X^ will 
thus be dependent, of lower order, or both. 

The importance of involutive symbols arises from the following theorem 
which provides a criterion for involution that can be tested in a finite num- 
ber of operations: 

Theorem 1 X is involutive if and only if Sym X is involutive and x[^ = X. 

10 Cartan-Kuranishi Completion 

The central theorem of the theory of involutive sytems of differential equations 
is the Cartan-Kuranishi theorem: 

Theorem 2 For every differential equation X there are two integers j, k such 

(k) 

that Xj is an involutive equation with the same solution space. 

The Cartan-Kuranishi completion algorithm is a straightforward application 
of the two previous theorems: 

input X 
repeat 

while SymZ is not involutive repeat X := X^ 

while X X^p repeat X := x[^ 
until Sym X is involutive and X = 
output involutive X 

Therefore, given a differential equation we may first complete it to an involutive 
system (if it is not already involutive) and then construct a power series solution 
order by order using the algorithm described earlier. 

11 Conclusion 

Although the algorithms described in this paper may be used to construct formal 
solutions to systems of partial differential equations, they suffer from several 
shortcomings in practice. Firstly, there is the problem of setting the values of the 
parametric derivatives. Typically, these must be calculated from the values of 
functions on submanifolds of the base space or from values of those functions at 
the points of a lattice within the base space. Secondly, many terms of the power 
series must be calculated to provide solutions of comparable accuracy to those 
produced by discretisation schemes, and the symbolic manipulations involved 
rapidly become computationally intensive as the order increases. Thirdly, there 
is the problem of the convergence of the series. A promising approach to the 
circumvention of these difficulties is the use of a hybrid method that first uses 
a discretisation scheme to calculate values at lattice points, uses these values to 
determine approximate values of the parametric derivatives and then constructs 
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power series about each of the points and smoothly joins them together using a 
functional interpolation scheme. This method is currently being implemented. 

Completion to involution and the construction of power series solutions are 
far from the only applications of the jet bundle formalism. As mentioned in 
the introduction, jet bundles also provide the natural setting for the analysis 
of the symmetry groups of systems of PDEs and of the variational symmetries 
of Lagrangian systems (which are linked by Noether's theorem to conservation 
laws). Symmetry analysis is also closely related to the construction of solutions 
possessing specified symmetries. Unfortunately, these fascinating and important 
subjects are beyond the scope of the current paper. 
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